Sound damping in ferrofluids: 
Magnetically enhanced compressional viscosity 
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The damping of sound waves in magnetized ferrofluids is investigated and shown to be considerably 
higher than in the non-magnetized case. This fact may be interpreted as a field-enhanced, effective 
compressional viscosity - in analogy to the ubiquitous field-enhanced shear viscosity that is known 
to be the reason for many unusual behavior of ferrofluids under shear. 
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I. INTRODUCTION 

Ferrofluids |Q are colloidal suspensions of mono- or 
subdomain ferrimagnetic nano-sized particles suspended 
in a carrier liquid. Under the influence of an exter- 
nal magnetic field the fluid behaves paramagnetically. 
Among the more remarkable flow phenomena of ferroflu- 
ids [Q, ^, ^ are the enhanced effective shear viscosity in 
a static magnetic field 0], or the viscosity decrease in 
response to an AC-field p[ 0, @] ■ Both are due to the so- 
called magneto-dissipative effect, which occurs when the 
experimental time scale compares to the magnetic relax- 
ation time. In those situations the actual magnetization 
M deviates significantly from its equilibrium value M cq . 
Then the increment M — M cq feeds back to the linear 
momentum balance, via the magneto-viscous stress ele- 
ment 
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M cq ) 



(1) 



leading to the appearance of an enhanced shear viscosity. 

When dealing with compressible flow situations such 
as sound, one needs to go beyond the approximation of 
incompressibility. If sound propagates through a magne- 
tized ferrofluid, density oscillations couple to the magne- 
tization, and one expects magneto-dissipation to become 
relevant as well. Sometimes the attenuation of sound 
is attributed to the elevated shear viscosity addressed 
above J8| , but this point of view disregards the fact that 
Eq ([!]) only contributes in shear flow geometries, and not 
in the compressional flows characteristic of sound: Tak- 
ing the divergence of the momentum balance (to derive 
an equation for V • v, the divergence of the velocity field) 
eliminates the contribution of ([!]), since ViVjAIIy = 0. 
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So if (|l|) were the only magneto-dissipative term, one 
must conclude that sound in magnetized ferrofluids does 
not experience any additional damping, but this is incor- 
rect: The recently derived ferrofluid dynamics S con- 
tains, in addition to the expression of Eq. (^), also a 
new, diagonal magneto-viscous stress element, which ac- 
counts for the additional energy loss of sound waves if 
the medium is magnetized. It is natural to interpret this 
fact as a magnetically enhanced compressional viscosity, 
in close analogy to the magnetically enhanced shear vis- 
cosity first observed by McTague |J] . Note that although 
this term was derived in || as a stringent result of energy 
and momentum conservation, it is not contained in the 
standard ferrofluid-dynamicsjl], ||, ||. 

While the well-known tensor element of (|l|) is non- 
vanishing only if the deviation (M — M Gq ) and the field 
H point in different directions, the new, diagonal stress 
element — writable as H • (M — M cq ) Sij in the case of 
linear constitutive relationship — remains non-vanishing 
even if both are parallel to each other. This is exactly 
the situation characteristic for the propagation of sound. 
Provided the sound frequency u> does not greatly ex- 
ceed the inverse magnetic relaxation r of the ferrofluid, 
a perceptible extra damping is predicted, several or- 
ders of magnitude larger than estimated by previous 
works fu], [jj. 

It is worth pointing out that the new diagonal stress 
element may of course be disregarded in incompressible 
flow situations. In these cases the pressure p{r,t) is de- 
termined by the condition V • v = 0. A diagonal stress 
element such as the above term then re-normalizes p, but 
leaves the velocity profile v(r, t) unchanged. Considera- 
tion of flow configurations such as those in sound, on the 
other hand, require the elimination of the incompress- 
ibility condition. This is done by adding the continuity 
equation for mass and adopting a viscous term propor- 
tional to V • v in the Navier-Stokes equation. In ordinary 
liquids, the pressure is now determined by the thermo- 
dynamic equilibrium relation p = p(p, T) as a function 
of density p and temperature T. But this is less simple 
in electrically polarized or magnetized liquids, where the 
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concept of pressure is rather ill-defined ||, |Tl|, ^2). The 
fact that the energy spent (or gained) by compressing a 
magnetized ferrofluid depends on the direction at which 
the force is applied implies that one needs to employ the 
more general concept of stress, and must especially be 
careful in handling the diagonal stress, as will be out- 
lined below. 

Several recent investigations on sound propagation 
in magnetizable fluids follow a more mesoscopic ap- 
proach, considering relative particle motions within clus- 
ters, aggregates or chains. Taketomi |l3| attributes the 
anisotropy of the sound attenuation coefficient to two 
types of motions performed by the ferrous colloidal par- 
ticles in the fluid, rotational and tranlational. Nahmad- 
Molinari et al. [jl4| investigated the propagation of 
collective modes through a magnetorheological slurry 
(with micron-sized grains) and observed two independent 
modes. The slower one, with a large amplitude, was con- 
sidered to be a compressional mode similar to what is 
found in porous fluid-saturated media. Later, Brand and 
Pleiner p"5| attributed this mode to a wave propagating 
along the particle chains. 

Others authors pursue a macroscopic, hydrody- 
namic treatment similar to the approach employed here. 
Parsons considered ferrofluids subject to strong mag- 
netic fields and took the vector of saturated magnetiza- 
tion as similar to the director in nematic liquid crystals. 
As a result, he found that magnetically induced rela- 
tive corrections to the sound velocity are small, around 
1CP 5 . As discussed by Henjes |nj, he worked with a 
purely mechanical pressure ignoring electromagnetic con- 
tributions in the diagonal stress. Using proper hydro- 
dynamics, Henjes pT| found that the magnetically in- 
duced corrections to the sound velocity are small, again 
of order 10 -5 . She also argued that since typical mag- 
netic relaxation times are of the order of r « 10 _6 s, 
magneto-dissipation is small for acoustic sound frequen- 
cies up to £j/(27r) = 20 kHz. This is correct, but it does 
not imply that magneto-dissipation can be entirely ig- 
nored, because all dissipation mechanisms derive from 
fast characteristic times, and the question is one of rela- 
tive weight. Summarizing, previous theoretical investiga- 
tions on sound propagation in ferrofluids do not account 
for magneto-dissipation. The present manuscript does, 
and the result is: In the hydrodynamic frequency regime, 
lot ~ 1, even moderate magnetic fields will induce extra 
damping of approximately 10%. 



II. THE STARTING EQUATIONS 

To quantify the damping in compressional flow sit- 
uations, propagation of sound waves through homoge- 
neously magnetized ferrofluids will be investigated. To 
streamline the consideration and to focus on the basic 
physics, we shall implement the following simplifications: 
(i) Only the leading order magnetic field effect 0(H 2 ) 
on the attenuation of sound is considered. This espe- 



cially implies the linear constitutive relation, M cq = %H. 
Moreover, the complication |l7j that sound in magnetized 
ferrofluids is generally accompanied by shear waves need 
not be considered. Although this coupling gives rise to 
rather surprising phenomena [T^ |, it contributes only at 
0(H 4 ) to the dispersion of sound (see below), (iii) We 
consider sound propagation and shear diffusion in the 
adiabatic limit. Adiabaticity means 5s = S(s/p) = 0, 
rather than 8T = as is the case in the isothermal limit. 
Adiabaticity is valid because in ferrofluids, shear diffu- 
sion and sound are usually fast processes on the time 
scale of heat conduction: The Prandtl number P, given 
by the quotient of characteristic thermal diffusion time 
over viscous diffusion time, or equivalently, by kinematic 
viscosity over heat diffusivity, P = v/k., is usually of 
the order of 10-100. (Depending on the ferrofluid, we 
have v w 10~ 6 - 10~ 3 to 2 /s, and k « 10~ 7 - 10^ 5 m 2 /s.) 
The same argument holds for solutal diffusion processes 
which are slower than shear diffusion and sound by a fac- 
tor P/L, where L is the Lewis number. For ferrofluids 
we typically have L = e>(10~ 4 ). 

Below it will be discussed in more details that the mag- 
netic susceptibility x> usually taken as a function of T 
and p, must then be considered as a function of entropy 
per unit mass s, in addition to p. Adiabaticity is a valid 
approximation here because in ferrofluids, shear diffusion 
and sound are usually fast processes on the time scale of 
heat conduction, (iv) Sound waves up to the MHz-range 
are weakly damped. The spatial decay length a -1 of the 
complex wave number k = lo/c + ia exceeds the wave 
number 2ttc/lo by many orders of magnitude 18 . Un- 



der those circumstances it is the custom to account for 
all damping mechanisms to linear order, (v) This paper 
focuses on sound attenuation. The tiny correction to the 
sound velocity is disregarded. An order of magnitude 
estimate for the magnetically induced correction yields 
Ac ~ (poxH 2 / p) 1 / 2 ■ Even at the highest magnetic field 
strength considered here, one gets Ac/c < 10~ 4 . 

The unperturbed state of the ferrofluid is given by a 
homogeneously magnetized ferrofluid at rest, with den- 
sity p and equilibrium magnetization M cq = ^H, where 
X is the magnetic susceptibility. To describe small am- 
plitude sound excitations we introduce deviations form 
this state Sp, v, <5H, (5B and <5M proportional to a plane 
wave with wave vector k. In particular the velocity field 
is taken as a longitudinal sound mode in the form 



V0C ]S e i(k.r+«t). 
k 



(2) 



The equations of motion governing the ferrofluid dynam- 
ics have recently been derived on the basis of the conser- 
vation laws and symmetries ||. The density field p{r,t) 
obeys as usual the continuity equation 



d t p + Vj{pVj 



0. 



The equation for the magnetization reads 



dt 



Mi - XiMiVi. 
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where d/dt = d t + v • V and fl = (V x v)/2 is the 
vorticity of the flow. The contributions proportional to 
Ai and A2 appear with the applied field breaking the 
isotropy of the system. These two terms reflect the fact 
that - in addition to the vorticity ft - compressional 
and elongational flow, denoted respectively as va = V • 
v and vfj — + VjVi — §<%i>fefc), contribute to 

the dynamics of M. Further terms associated with the 
uniaxial symmetry (see the terms proportional to A3 and 
A4 in Eq. (13) of Ref.j|]) have been omitted on the left 
hand side of Eq.(^), as they are of higher order in the 
magnetic field. The increment 



coefficients Ai and A2 (not yet measured) should be neg- 
ligibly small. As outlined in Sec. fl this term arrives 
cogently during the derivation of the stress tensor and 
accounts for magneto-dissipative processes if h is paral- 
lel to the equilibrium magnetization M ei? . This is crucial 
for situations where M and H oscillate co-linearly but 
with a temporal phase lag. (Recall that the customary 
magneto-dissipative term given by Eq. (Q) drops out if 
M and H are parallel to each other.) 



III. RESULTS 
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with B cq = /xoM(l + X)/Xi accounts for magneto- 
dissipative relaxation, on the time scale given by r. 

The magnetic field variables H, B, M are defined in 
Si-units as usual, with po the vacuum permeability. For 
the evolution of the magnetic fields we adopt the static 
Maxwell equations V-B = VxH = 0. With the plane 
wave behavior similar to Eq. (g) the fluctuations are re- 
lated in the following manner 



£H = -<SM|i; 5B = f i SM ± . 



(0) 



Here the indices || and _L refer to the respective directions 
relative to the propagation direction k of the wave. 

The balance equation for the linear momentum reads 
dtpvi + VjHj = 0, with the stress tensor |j| 

Hy = [-u + sT + up + H • B (7) 
-(Ai-^h-M^-n^-ff^ 

-^-(Mihj + Mjhi) + ^(hiMj - hjMi). 

To avoid misunderstandings of what is meant by the 
" pressure at non-zero magnetic field strength" , the di- 
agonal element is written in terms of the density of total 
energy u, the entropy density s, and the chemical poten- 
tial p. The viscous stresses IIy a = 2r\\ + 772 SijVkk are 
taken as usual with the shear viscosity 771 and the volume 
viscosity 772. The terms proportional to A12 are counter 
terms to those of Eq. (^), they are constrained by the 
Onsager symmetry relations. 

To make contact to previous formulations of the stress 
tensor |Q || we have to switch for a moment to T rather 
than s as an independent variable. Then the square 
bracket in Eq. (Q) can be recast in terms of the thermo- 
dynamic relation for the pressure at zero magnetic field 
Po{p,T) 

Po+Mo^-(Ax-iA 2 + l)h.M + ^(l-^ (8) 
2 3 2 X X dp 

where x — xiP^T). Note that magneto-dissipation, pro- 
portional to h • M, remains finite even if the transport 



A. Dispersion of isothermal sound waves 

We now return to the adiabatic formulation, where x = 
x{p,s). Using Eq.(|^) and the longitudinal plane wave 
velocity field (Q) the magnetization fluctuation <5M = 
<5Mii + <5Mj_ is related to the density variations by 



<5M|| = M 



[ {p/x)dx/dp-iuT{\ 1 
\ 1 + X + iur 
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The real and imaginary part of <5M is associated with 
magnetically induced corrections to the sound velo city 
and attenuation, respectively. Substituting Eqs.(p||ic|) 
into the divergence of the momentum balance yields the 
following complex dispersion relation for sound waves in 
magnetized ferrofluids 



C s 



2pci 



(11) 



where c 2 s = dpo(p, s)/dp is the square of the zero-field 
adiabatic sound velocity. Recall that - according to ap- 
proximation (v) - magnetic corrections of c s are disre- 
garded. The increment r\ m is given by 



Vm = PotxH 
where 



k± sin 



K 



(l + x) 2 + (™) 2 l + (m,) 2 



-^ + (1 + x)(Ai + |a 2 ) 
X op 3 



X op 



(12) 



(13) 



(14) 



i] m can be interpreted as a "magnetic extra viscosity". 
The expression for rj m is clearly anisotropic, and 9 de- 
notes the angle between the applied magnetic field H 
and the direction k of propagation. Note also that rj m is 
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frequency dependent, being maximal at tlo — > and van- 
ishing in the high frequency limit tlo ^> 1. Eqs.([l3 14) 
can be simplified if the magnetic susceptibility is pro- 
portional to the density thus {pjx)dx/dp w 1. For a 
rough estimate, take the data Il9[ for a high viscosity 
hydrocarbon based ferrofluid (APG 933, Ferrofiuidics) : 
X — 1.1, Tji ~ 0.5Pas, and r ~ 0.55ms, in addition to 
Ai = A2 = r/2 = (for lack of better information). Then 
a lOOHz-sound wave propagating in an applied magnetic 
field of say H = 10 4 A/m (this is a field strength at which 
most fcrrofluids still obey linear constitutive relations) 
experiences a magneto-viscous extra damping of ~ 7% 
at the parallel orientation k |j H (i.e. 6 = 0) and almost 
9% at the transverse setup (0 = 90°). The same estimate 
applies to a ferrofluid of similar microscopic make-up, but 
at a viscosity of r]\ ~ 5 mPas and a frequency of 10kHz. 
(Here we assume r cx rji, valid for Brownian particles, 
i.e. when the particle's magnetic moment is fixed to 
the crystallographic orientation). Damping increments 
of this size should be detectable in a careful sound wave 
experiment. Moreover, by scanning the ^-dependence of 
r] m , it should be possible to obtain information on the 
transport coefficients Ai and A2. 



B. Adiabatic versus isothermal susceptibility 

Knowing the dependence of the magnetic susceptibil- 
ity as a function of density and temperature, x{p^), 
the derivative p^(p,s) in the above equations can be 



expressed as follows 
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(15) 



involving both the magneto-strictive and the magneto- 
caloric contributions. Here a v = ~(l/p)dp(p,T)/dT de- 
notes the thermal expansion coefficient, ct the isothermal 
sound velocity, and C v = Tds(T, p) / dT the specific heat 
at constant volume, everyone of them evaluated at zero 
magnetic field. Eq. ([if]) indicates that adiabatic sound 
waves involve both magneto-strictive as well as magneto- 
caloric contributions. Here a v = —{\/p)dp(p,T)/dT 
denotes the thermal expansion coefficient and C v = 
Tds(T, p)/dT the specific heat at constant volume, each 
of them evaluated at zero magnetic field. For a typi- 
cal olefine-based carrier liquid the dimcnsionlcss factor 
c^a v /C v can be estimated by 0.3. 



C. Enhanced compressional viscosity 



In order to classify the viscosity increment rj n 



[Eq. (12)] as an field dependent offset to either the shear 
viscosity 771 or the volume viscosity 772 we evaluate the en- 
tropy production in the present setup. Following Ref. |J 
the total entropy production is given by 



Computing the magneto-viscous surplus [last term in 
(|l6|)1 up to first order in lot yields 



l2 rj2 

h = p. rxH 

POT 



(1 + x) 2 
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The formal similarity of (|l^ ) with the second term of 
( filf ) suggests that the "magnetic extra viscosity" r\ m ac- 
cording to Eq. (|lj) is to be interpreted as an enhanced 
compressional viscosity Ai]2(H). 



D. Comparison with experiments 

The experimental material on sound propagation in 
fcrrofluids is rather scarce. The early measurements of 
Chung and Isler ]|, |2^] on a water-based ferrofluid seem 
to be the only available systematic study of the veloc- 
ity and attenuation of sound (note, however, that Sku- 
miel's ||l| later investigations reveal a strong dependence 
of the sound velocity on the type of the carrier liquid, 
i.e whether it is aqueous or organic). The experiments 
of Refs. H |T| were carried out with 2.25 MHz ultra- 
sound, employing pulse-echo and continuous wave meth- 
ods. The experimental data cover a wide magnetic field 
range from up to 2500Gauss = 2 x 10 5 A/m. Within 
the weak field subrange, where linear constitutive rela- 
tions hold for most ferrofluids, H < 10 4 A/m= 125Gauss 
say, the damping increment a was found to increase by 
1.8dB (~ 20%) at = but to decrease (anomalous 
sound attenuation) by almost 3.5dB (~ 50%) at 6 = 90°. 
Unfortunately no information was given whether demag- 
netization effects due to the cylindrical probe geometry 
had been taken into account here. However, we point 
out that the observed anomalous -ff-dependence does 
not even qualitatively comply with the present theory. 
Neither the observed history dependence of the exper- 
imental data (also detected by Gotoh et al. f^l ) can 
be explained by the present approach. The latter pecu- 
liarities suggests that the small ultrasound wavelengths 
couple to microscopic inhomogeneities associated with 
particle chains or clusters in the ferrofluid suspension. 
Anisotropies, field dependencies and anomalies recorded 
by ultrasound experiments therefore seem to depend on 
mechanisms which are rather different from those covered 
by the present hydrodynamic analysis. Owing to the lack 
of other pertinent experimental data let us nevertheless 
try a quantitative comparison with the measurements of 
Isler and Chung ||. For a rough estimate we take the 
viscosity of their aqueous ferrofluid by 771 = 10~ 3 Pas 
and the susceptibility by \ ~ 1 (specifications are not 
given). Due to the high ultra-sound drive frequency the 
experiment was operated in the limit lot 3> 1 and thus 
the expected extra damping is quite small. Even if the 
magnetic relaxation time r is estimated to be as small 
as 10~ 6 s, the prediction of Eq.(p"2|) at 6 = is two or- 
ders of magnitude smaller than the empiric value. We 
therefore conclude that for a reliable quantitative check 
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of the present theory, experiments at acoustic frequencies 
(where lot ~ 1) would be more suitable. 



IV. DISCUSSION 



Assuming that a sound emitter produces plane density 
waves Sp(t) within a magnetized ferrofluid, the right hand 
side of Eq. (jl|) can be recast as 



The present analysis deals with the attenuation of 
sound in fcrrofluids, which are exposed to a weak ho- 
mogeneous magnetic field in any direction relative to the 
propagation. This has been accomplished by investigat- 
ing the linear dispersion of a pure longitudinal veloc- 
ity excitation. Recently, it has been pointed out fTjj ] 
that density excitations (sound) and transverse velocity 
fluctuations (shear waves) in magnetized ferrofluids do 
not evolve separately as is the case at H = 0: At fi- 
nite H, sound waves may produce shear excitations and 
vice versa. Clearly, if shear waves accompany sound this 
opens a new attenuation mechanism, which cannot be ig- 
nored. In the remainder of this section we shall argue why 
this cross-coupling remains without consequences for the 
present analysis: By taking the curl of the momentum 
balance one arrives at 

pdtQ. - r;V 2 = -VxVx(hx M). (18) 



V x V x (h x M) 



rpo d\ M x x M || 
X dp 



X 



V 2 d t 6p 



A 9 ) 

thus acting as a magneto-dissipative source of vorticity 
O. If the applied magnetic field is weak we have Q — 
0(H 2 ). Via the term $7 x M in Eq.(j|) this sound-made 
vorticity induces a third order correction in M, which - 
at the considered accuracy level 0{H 2 ) - does not affect 
the sound dispersion. Note however, that a proper study 
of sound damping in strongly magnetized ferrofluids (as 
for instance undertaken in Ref. |10fl ) must not ignore the 
complication arising from the magneto-dissipative cross- 
coupling between compressional and shear excitations. 
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